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Abstract 

The entropy of a polymer confined in a curved surface and the elastic free energy of a membrane 
consisting of polymers are obtained by scaling analysis. It is found that the elastic free energy 
of the membrane has the form of the in-plane strain energy plus Helfrich's curvature energy [Z. 
Naturforsch. C 28, 693 (1973)]. The elastic constants in the free energy are obtained by discussing 
two simplified models: one is the polymer membrane without in-plane strains and asymmetry 
between its two sides, which is the counterpart of quantum mechanics in curved surface [Jensen 
and Koppe, Ann. Phys. 63, 586 (1971)]; another is the planar rubber membrane with homogeneous 
in-plane strains. The equations to describe equilibrium shape and in-plane strains of the polymer 
vesicles by osmotic pressure are derived by taking the first order variation of the total free energy 
containing the elastic free energy, the surface tension energy and the term induced by osmotic 
pressure. The critical pressure, above which spherical polymer vesicle will lose its stability, is 
obtained by taking the second order variation of the total free energy. It is found that the in- 
plane mode also plays important role in the critical pressure because it couples with the out-of- 
plane mode. Theoretical results reveal that polymer vesicles possess the mechanical properties 
intermediate between fluid membranes and solid shells. 

PACS numbers: 61.41. +c 



*Email address: oy@itp.ac.cn 



2 



I. INTRODUCTION 



Thin thickness structures exist widely in Nature. Many things such as eggs, snails, air- 
planes and so on in our daily life are covered with solid shells that play protective roles. 
In the realm that we cannot see with naked eyes, virus usually have protein shells, and 
eukaryotic cells are enclosed by cell membranes that consist of lipids, proteins and carbohy- 
drates etc. A lipid molecule has a polar hydrophilic head group and one or two hydrophobic 
hydrocarbon tails. When a quantity of lipid molecules disperse in water, they will assem- 
ble themselves into a lipid bilayer in which the hydrophilic heads shield the hydrophobic 
tails from the water surroundings because of the hydrophobic forces. Solid shells and lipid 
bilayers are, respectively, in the categories of hard condensed matter and soft one. Their 
mechanical properties have attracted much attention for a long time □ DO 

The significant difference between solid shells and lipid bilayers is that the former can 
endure the in-plane shear stress but the latter cannot. Due to this difference, solid shells 
and lipid bilayers have different forms of deformation energy. Under the assumption of 



homogenous and isotropic bulk materials and in the 



energy per unit area of a solid shell is expressed as [l]: 

£ sh = | [{2Hf - 2(1 - u)K] 

+ ^r^y [( 2 J? - 2(1 - m , (i) 

where D = (l/12)Yh 3 /(l — v 2 ) and C = Yh are bending rigidity and in-plane stiffness of 
the shell. Y, v and h are, respectively, the Young's modulus, the Poisson ratio and the 
thickness of the shell. 2 J and Q are the trace and determinant of the in-plane strain tensor, 
respectively. For a spherical solid shell with radius R, the critical osmotic pressure (i.e., the 
pressure difference between the outer surface and inner one of the shell, above which the 
shell loses its stability) is 

_ 2Yh 2 

PCS ~ V3(l-^ 2 ' ( } 

In 1973, Helfrich [3j recognized that the lipid bilayer was just like liquid crystal in smectic 

A phase at room temperature. Based on the elastic theory of liquid crystal 7(, he proposed 

the curvature energy per unit area of the bilayer 

£ lb = (k c /2)(2H + c ) 2 + kK, (3) 



imit of thin thickness, the elastic free 



where k c and k are elastic constants; and H, K and cq are mean curvature, Guassian 
curvature and spontaneous curvature of the lipid bilayer, respectively. For phospholipid 
bilayers at room temperature T, the persistence length is usually much larger than the 
size of the membranes and the effect of shape fluctuations is negligible 4] because of k c « 
10~ 19 J^> ksT where ks is the Boltzmann factor. The free energy of a closed bilayer 
under the osmotic pressure p is written as = §{£u> + X)dA + p J dV, where dA is the 
area element and V the volume enclosed by the closed bilayer. A is the surface tension of 
the bilayer. The first order variation of T\h gives the shape equation of closed bilayer 

p - 2XH + k c V 2 (2H) + k c (2H + c )(2H 2 - c H - 2K) = 0. (4) 

For a spherical lipid bilayer with radius R, the critical osmotic pressure for stability is [10]: 

2k c {6 - c R) _ 
Pd = ^3 • (5) 

It follows that pd ~ k c / R? because the typical value of c R is about 1. Therefore, lipid 
bilayer is, indeed, much softer than solid shell. 

Are there membranes intermediate in the state between Helfrich's fluid lipid bilayers and 
classical solid shells? The polymer vesicles discussed below may be as a example. In the 
last decade, Decher invented the layer-by-layer assembling technique Following this 

technique, Caruso et al. made spherical polyelectrolyte capsules by the step-wise adsorption 
of p olvelectrolytes onto charged colloidal templates and then decomposition of the templates 
12l 1 1 31 ] . The capsules composed of about 10 layers of alternating polystyrene sulfonate and 
polyallylamine hydrochloride. The thickness h of capsules was about tens of nanometers 
which was remarkably less than their radii R (several micrometers). Gao et al. jlj] found 
that the spherical polyelectrolyte capsule lost its stability and changed its shape abruptly 
above some threshold of osmotic pressure p c which was proportional to R~ 2 and h 2 . In 
their experiment, the thickness dependence of p c might not be exact because the polyelec- 
trolyte capsule with more than 10 layers was chemically instable as they claimed. Thev 
also explained their results through the stability theory of classical elastic solid shells |l|, |2| . 
But it is well known that the classical theory is based on the assumption of homogenous 
and isotropic bulk materials which entirely ignores the characteristic of the polyelectrolyte 
capsule consisting of many polymers. If considering the polymer structures of spherical 
polyelectrolyte capsule, can we still derive p c ~ i? -2 ? 



In this paper, we will answer above questions. To do that, we derive the entropy of a 
polymer confined in a curved surface and the elastic free energy of a membrane consisting 
of polymers by scaling analysis. It is found that the elastic free energy of the polymer 
membrane has the form of the in-plane strain energy plus Helfrich's curvature energy. The 
elastic constants in the free energy are obtained by discussing two simplified models: one 
is the polymer membrane without in-plane strains and asymmetry between its two sides, 
which is the counterpart of quantum mechanics in curved surface [l5|; another is the planar 
rubber membrane with homogeneous in-plane strains. The equations to describe equilibrium 
shape and in-plane strains of polymer vesicles by osmotic pressure are derived by taking the 
first order variation of the total free energy containing the elastic free energy, the surface 
tension energy and the term induced by osmotic pressure. The critical pressure, above which 
a spherical polymer vesicle will lose its stability, is obtained by taking the second order 
variation of the total free energy. It is found that the in-plane mode also plays important 
role in the critical pressure because it couples with the out-of-plane mode. These theoretical 
results reveal that polymer vesicles possess the mechanical properties being intermediate 
between Helfrich's fluid membranes and classical solid shells. 

The following contents of this paper are organized as below: In Sec. |Hj we derive the free 
energy of polymer membrane by using the scaling concepts [16i]. In Sec. IHIl we obtain the 
shape and in-plane strain equations of closed polymer vesicles by using the surface variation 



theory developed in Ref. JjJ, In Sec. IIV[ we discuss the mechanical stability of the 
spherical polymer vesicle by taking the second order variation of the free energy. In SecfV] 
we give a brief summary and discussion. 



II. THE FREE ENERGY OF POLYMER MEMBRANE 

The polymer membrane concerned in this paper is one or a few thin layers consisting of 
cross-linking polymer structure like rubber [l9| at molecular levels. It can be represented as 
a mathematical surface with curvature and strains. It is hard to derive its free energy in a 
strict way. But we can drive it by using scaling concepts in polymer physics proposed by de 
Gennes [16]. In the following contents, we take de Gennes' convention: the entropy S is a 
dimensionless quantity and Boltzmann factor ks is implicated in temperature T. 



5 



A. The free energy of polymer membrane 

If we take the Gaussian chain model j^J, the root of mean square end-to-end distance of 
a polymer is Rq ~ y/Nb®, where b$ is segment length of the polymer and N is the number 
of segments. Assume that the principal radii of the surface are much larger than Rq. If 
the in-plane strain tensor of the surface is denoted by e which is assumed to be a small 
quantity, the entropy of the polymer confined the surface must be the function of 2HR , 
KRq, 2 J and Q because it is a dimensionless invariant quantity under the transformation of 
coordinates, where H, K, J = tre and Q = tie are the mean curvature of the surface, the 
Gaussian curvature of the surface, the trace of strain tensor and the determinant of strain 
tensor, respectively. Thus we can expand it as 

S ~ A^HRq) + A 2 (2HR ) 2 + A 3 KR 2 + B 2 {2J) 2 + B 3 Q (6) 

up to the second order terms, where A±, A2, A3, B 2 , B 3 are constants. In this expression, an 
unimportant constant term is neglected. Moreover, one must notice that generally we have 
K ~ (2H) 2 and Q ~ (2J) 2 . Additionally, there is no first order term of 2J in the expression 
of the entropy because we expect that — e plays the same role as e in the entropy. It is useful 
to write the entropy in another equivalent form 

S ~ A 2 R 2 (2H + c ) 2 + A 3 R 2 K + B 2 (2J) 2 + B 3 Q, (7) 

where c = Ai/{2A 2 ) is a constant, so called spontaneous curvature, which is expected to 
satisfy |co-Ro| ^ 1- 111 fact, cq vanishes if there is no asymmetric factor between two sides of 
the surface because H turn into —H if we change the normal direction of the surface. 

Assume h to be the membrane thickness and M the number of polymers per volume. 
Additionally, we neglect the the entanglement of polymers. Consequently, the free energy 
per unit area of a membrane consisting of polymers has the following form 

S pm = —(Mh)TS 

= ^[(2J) 2 - nQ] + |(2tf + c ) 2 - kK, (8) 

where k d = —2B 2 MhT, /i = -B 3 /B 2 , k c = -2MhTA 2 R 2 , k = MhTA 3 R 2 . 

Obviously, Eq.( | 8* ]) is degenerated to Helfrich's curvature energy of fluid membranes if 
k d = 0, and to the elastic energy of classical solid shells if c = 0, kd = C/(l — u 2 ), k c = D, 
H = 2(l-v) and k = D(l-v). 
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B. The elastic constants kd, fi and k c 

kd, /i, k c and k are unknown universal constants independent of the detailed shape and 
the small in-plane strains of the polymer membrane. If only discussing the closed polymer 
vesicles in this paper, we need not to know k because the integral of kK is an unimportant 
constant. To determine kd, p and k c , we will discuss two ideal cases: One is the cylinder 
polymer membrane without any strain and asymmetry between its two sides; another is the 
planar membrane with the homogenous in-plane strains. 

In the former case, we denote p the radius of the cylinder. On the one hand, (jH)) is 
simplified as 

k c 

On the other hand, we know there is a 1-1 correspondence relation between polymer statis- 
tics and quantum path integral method 20] as shown in Table UJ In 1971, Jensen and 
Koppe dealt with the quantum mechanics of a particle constrained in a curved surface and 
obtained a nontrivial conclusion [l^ : the constraint would induce an effective potential 
V e f = — ^[(2H) 2 — AK] in Shrodinger equation, where m is the particle mass. In terms of 
the correspondence rules in TableHJ there will be an effective potential U e f = -^[(2H) 2 — AK] 

Tb 2 

for a polymer confined in curved surface. Especially, U e f = -^z for the cylinder with radius 
p. We must pay more attention to the fact that there is a minus symbol in the potential term 
when we use the correspondence rules. In fact, Yaman et al. overlooked this fact in recent 
literature [2]]. But this flaw can not diminish the value of their pioneer work in the study 
of polymer confined in curved surface. Thus their results can be safely transplanted only if 
we change the sign. Consequently, we obtain the free energy of the cylindrical membrane 
consisting of Mh polymers per unit 

MhTNbl _ MhTR 2 
pm 24p 2 24p 2 1 ' 

if neglecting the entanglement between polymers. Comparing Eq.fjHJ) with Eq.()10|). we obtain 
k c = MhR 2 T/12. 

In the latter case, H, K and Cq are vanishing for planar membrane with symmetry between 
its two sides. On the one hand, (JEJ) is simplified as 

S vm = ^[{2Jf-pQ}. (11) 
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For homogenous stain e, we can express it by its components en, 622 and 612 = 621 = in 
some orthonormal coordinate system so that 2 J = en + 622 and Q = 611622- 

On the other hand, we notice that there might be cross-linking joints between polymers in 
the membrane. This character suggests that the membrane should have the elastic properties 
of rubber materials. In terms of the elasticity theory of rubber h| , the deformation energy 
of a planar rubber per area can be expressed as f r = (MhT/2) [Af + A| + 1/ (\ 2 \ 2 ) — 3] , where 
Ai = 1 + en and A2 = 1 + 622 are extensions. For small strains, it is expanded to the lowest 
order terms as 

S pm = 2MhT(e 2 u + 611622 + ej 2 ) = 2MhT[(2J) - Q}. (12) 

Thus we can obtain = AMhT and /i = 1 by comparing Eq. (jllj) with Eq. (jl2|) . 

Therefore, we obtain the free energy of a closed polymer vesicle under osmotic pressure 

p: 

T= j{k d /2)[{2J) 2 -Q]dA + j[{k c /2)(2H + c ) 2 + y\dA + p J dV, (13) 

where fi, A and V are surface tension, surface area and volume enclosed by the vesicle, 
respectively. In this expression, the term related to the Gaussian curvature disappears 
because its integration § KdA is an unimportant constant so that it is omitted. It is easy 
to see that Eq. (fT3^) is degenerated to the free energy of closed lipid bilayer for ka = 0, and 
to the free energy of solid shell with v = 1/2 if c = 0, = C/(l — v 2 ) and k c = D. 



III. THE SHAPE AND IN-PLANE STRAIN EQUATIONS OF CLOSED POLY- 
MER VESICLES 

In this section, we will give the shape and in-plane strain equations of closed polymer vesi- 
cles from the first order variation of free energy (|13|) . The method has been fully developed 
in Ref . , and the key elements and notations are shown in the Appendix. 

If a point ro in a surface undergoing a displacement u to arrive at point r, we have 
du = dr — dr and naturally Sidu = 5idr (i = 1, 2, 3). 

If denote dr = uiei + co> 2 e 2 and du = UiUi + U2CU2 with |Ui| 1, |U 2 | 1, we can 
define the in-plane strains [22]: 



£11 



du ■ ei 
IdrJ' 



« Ui ■ ei, (14) 

oj 2 =0 



8 



^22 



du ■ e 2 

\dr \ 



U 2 -e 2 , (15) 



wi=0 



1 

£12 



du ■ e 2 \ Z' du ■ ei 



CJX"n / n V ^I"n / 

11 ' aj2=0 \ I u l / ui=0 



~(Vi-e 2 + U 2 -e 1 ). (16) 



2 

Using <5j(iu = <5j(ir and the definitions of strains (|14 p — we can obtain the leading 



terms of variational relations: 

SiSiiUi A u 2 = SiUi A u 2 , (17) 

5ie 12 uJiAuj2 = - [oji A SiUx + SiU 2 A u 2 ] , (18) 

5i£ 2 2^i A oo 2 = uji A 5jCJ 2 . (19) 

From Eas. (IXT2|) - (|X2T| and (fT7 jl -(P |l . we have: 

XT i"; r J/on, £lld^2 - ^12^1 . d(ei2^1 +£22^2) ln /onN 

= <p k d [-d(2J) Au 2 + -i '-]n u (20) 



6 2 T = j> k d [d(2J) A Wl - £A ~ &2 A - d ( g n"i + g ^) ]n a> (21) 



= j [k c {2H + c )(2H 2 - c H - 2K) + k c W\2H) 

+p - 2H(n + k d J) - -| (aen + 2te 12 + ce 22 )]n 3 dA. (22) 
Thus the Euler-Lagrange equations corresponding to the functional (|T3J) are 

k d [-d(2J) Aw 2 - ^{suduj 2 - e 12 dui) + -d(ei 2 wi + £22^2)] = 0, (23) 

k d [d(2J) A a>i — -(ei 2 G?o? 2 - £ 2 2<M) - ^(eii^i + £i 2 u; 2 )] = 0, (24) 
p - 2#(/x + fc d J) + A; c (2# + c )(2F 2 - c # - 2tf) + k c V 2 (2H) 

-k( a£u + 2 be 12 + ce 22 ) = 0. (25) 

Eqs. f|23jl and (|24*|) are called the in-plane strain equations because they describe the in-plane 
strains of polymer vesicles under the pressure p. Eq.(J2SJ is called the shape equation because 
it describes the equilibrium shape of polymer vesicles under the pressure p. 

Obviously, if k d = 0, then Eqs. (|23p and (J2jj) are two identities while Eq.(|'25p is degen- 
erated into shape equation (JU) of closed lipid bilayers. Generally speaking, it is difficult to 
find the analytical solutions to Eqs. (J23j) - (j2l)j) . But it is easy to verify that en = e 22 = e (a 
constant), t\ 2 = satisfy Eqs.()23 |) - (j25|) for a spherical vesicle with radius R if the following 
equation is valid: 

pR 2 + (2/i + 3k d e)R + k c c {c R - 2) = 0. (26) 
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IV. MECHANICAL STABILITY OF SPHERICAL POLYMER VESICLES 



Now we will calculate the second order variation of functional (JT3j) and discuss the me- 
chanical stability of a spherical polymer vesicle. 

In Ref. [uj], only the term b\T related to the out-plane mode {fi 3 } is calculated. Here we 
also consider the contribution of in-plane mode {fii,!^}. Due to the notation of exterior 
differential d and Hodge star *, Qi and can be expressed as 

fii^i + Q 2 u 2 = dQ + *dx (27) 



by two scalar potential functions il and x f° r 2D manifold |23|. Using Eqs. ()A12|) - ()A21|) and 
(fT7|) - (Pj) . we can calculate 5\ b\T, 6\? \ S^F, Ms-? 7 , and S^J 7 from Eqs.fl2DJ-©2J and 
(f2T)j) for spherical polymer membrane. Eventually, we arrive at 



where 



5 2 F = b\T + b\T + b\T + IbxhT + Ibxh? + 25 2 5 Z F = G x + G 2 , (28) 

Gi = (j) tt 2 3 {3k d /R 2 + (2k c c /R 3 ) +p/R}dA 

+ j tt 3 V 2 tt 3 {k c c /R + 2k c /R 2 + P R/2}dA + j k c (\7 2 tt 3 ) 2 dA 

j tt 3 V 2 ttdA + k d j ( V 2 fi) 2dA+ ^j VV 2 ttdA, (29) 



3k d 
+— 
R 



G 2 = ^ <f(V 2 X fdA + x V X dA. (30) 

If we take k = k c /2, K = 3/^/2, fi = k^/2, w = Q 3 and \& = Q in equations (6) and (7) of 
Zhang et a/.'s paper j^j], then G\ and G 2 correspond to Fi[io ) \E f ] and F 2 [x] in that paper 
under the conditions of p = and CqR = 2. Obviously, there is no coupling between modes 
{x} and {f2,f2 3 }; but there is coupling between in-plane mode {fl} and out-of-plane mode 
{f^}. We will show that in-plane modes have quantitive effect on the stability of the cell 
membrane although they can not qualitatively modify the results of Ref. [18]. 

Because G 2 is obviously positive definite, we merely need to discuss G\. fl 3 and Q 
in the expression of G\ can be expanded by spherical harmonic functions 25| as = 



Yl™oY^m= l -l a lrn Y lm{Q,<l>) an d ^ = EtoS=-|Wm(^) witn a lm = (-l) m %-m and 

b* m = (-l) m b^_ m . It follows that 

00 I 

G x = ^^2|a im | 2 {3A; (i +[/(/ + l)-2][/(/ + l)A ;c / J R 2 -A; c co/ J R-j9i?/2]} 



1=0 m=0 
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3k n 



~ + l )( a *lm h l m + ai m b* lm ) 



R 

1=0 m=0 
oo / ; 

EE 

1=0 m=0 



+ LLi[ 2|2 ( | + 1 ) 2 - | ( | + 1 )]l^l 2 ( 31 ) 



We find that if p < Pl = m +§_ 1]R + M^£gg] (/ = 2, 3, • • •), then d is positive 
definite, i.e., the vesicle is stable. We must take the minimum of p\ to obtain the critical 



pressure: 



+ ^c^-coKj < M^coKj ; (3 ^ 2 < 121A . c 



p c = minfe} = { L1H w K ° ' v " . (32) 

2VpK + ^(i_ 2c0J R), (3^ 2 > 121fc c ) 

But if we do not consider the in-plane mode {^}, we will obtain the critical pressure [rsl ]: 



J» + < 2fce(1 ^r ofl) , (s^ 2 < i6fc e ) 

PC I ^p^ + |f(2-coi?), (3^ 2 >16fc c ) 

Comparing Eq. (J32|) with (J33|) . we find that in-plane mode have quantitive effect on the 
stability of the polymer vesicles although they can not qualitatively modify the result without 
considering it. 

Now we test the validity of Eq.(j32J) by considering two special cases. The first case, 
kd = 0, corresponds to lipid bilayer. In this case, Eq. (jH2^1 is exactly reduced to Eq.(jHJ), the 
critical pressure for spherical lipid bilayer. The second case, Co = 0, kd = Yh/(1 — z/ 2 ) 
and k c = Yh 3 /[12(1 — u 2 )] with v = 1/2, corresponds to the solid shell with Young's 
modulus Y, Poisson ratio u, shell thickness h. Under the condition of h <C R, Eq. (}3*2*|) gives 
p c = (4/3)Yh 2 / R 2 that is exactly the result of Eq.(J2J) with v — 1/2. Thus we are sure of the 
validity of Eq. (j32jl . 

Now we turn to the polymer vesicle consisting of polyelectrolytes. In the experiment by 
Gao et al. Q], its thickness is h ~ 20nm which is much smaller than its radius R ~ 2/im. 
The segment length is taken as 3 times of carbon-carbon bond length, i.e., bo ~ 4. 2 A. 
The number of segments per polymer are about N ~ 200 due to the molecule weight 
70000. Thus Rq ~ VNb ~ 60A, which is less less than R. Considering kd = AMhT and 
k c = MhR 2 T/12, we arrive at 3k d R 2 / (I21k c ) = (12R) 2 / (1LR ) 2 > 1, i.e., 3k d R 2 > 121k c . 
Under this condition, (|3*2|) is reduced to 

2MTR h 

Pc = ^ • (34) 
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Above equation can explain the experimental result p c ~ R~ 2 obtained by Gao et al. But 
their result p c ~ h 2 is inconsistent with our theoretical result p c ~ /i. The main reason for 
this discrepancy is that the membrane consisting of polyelectrolytes is chemically instable 
if its layer number is more than 10, which will disturb to test the exact relation between 
critical pressure p c and thickness h. 



V. CONCLUSION 



In above discussion, we briefly introduce the polymer statistics and the correspondence 
principle (shown in Table HJ) between it and path integral method in quantum mechanics. 
We derive the entropy of a polymer confined in a curved surface and the elastic energy of a 
membrane consisting of polymers by scaling analysis. It is found that the elastic energy of 
the polymer membrane has the form of the in-plane strain energy plus Helfrich's curvature 
energy as shown in Eq.(jHJ). The elastic constants kd, k c , \i in the free energy are obtained by 
discussing two simplified models: one is the polymer membrane without in-plane strains and 
asymmetry between its two sides, which is the counterpart of quantum mechanics in curved 
surface; another is the planar rubber membrane with homogeneous in-plane strains. The 
equations to describe equilibrium shape and in-plane strains of polymer vesicles by osmotic 
pressure are derived by taking the first order variation of the total free energy (|13j) containing 
the elastic free energy, the surface tension energy and the term induced by osmotic pressure. 
The critical pressure (|32jh above which the spherical polymer vesicle will lose its stability, is 
obtained by taking the second order variation of the total free energy ()13|) . It is found that 
the in-plane mode {Q} also plays important role in the critical pressure because it couples 
with the out-of-plane mode 

We estimate that p c = 2 M^oh through the experiment by Gao et al. This result is 
qualitatively intermediate between p c i = 2fce ( 6 ~ cofl ) f or lipid bilayer and p cs = , 2Yh ' 2 — for 

n y 3(1— v 2 )R 2 

solid shell. Therefore polymer vesicles possess the mechanical properties being intermediate 
between Helfrich's fluid membranes and classical solid shells. But is it reasonable to use 
the present theory to polyelectrolyte membranes? We discuss two points: (i) The derivation 
of the bending rigidity of the polymer membrane uses results for the quantum mechanics 
of a particle constraint in a curved surface. The results can only apply to polymers that 
are much more constraint perpendicular to the membrane than their lateral size. It seems 
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that the polymer vesicles in the experiment by Gao et al. do not satisfy this condition 
since Rq ~ 6 nm and ~ 20 nm. But in fact, the vesicles contain 10 layers and the 
thickness of each layer is about 2 nm that is much small than R and the total length 
of a polymer. That is, the strong constraint in the normal direction of the membrane is 
satisfied, (ii) The Gaussian chain model is used in the present work. However, the polymers 
in polyelectrolyte membranes may be nonideal. This is indeed a difficulty But the present 
theory is focused on the small deformations of polymer vesicles, the model of ideal polymer 
should give approximate results due to the lessons in classical theory of rubber elasticity 
[lflj ]. Additionally, we indeed obtain the relation p c ~ R~ 2 observed by the experiment. 

It is a nontrivial thing that we analysis the mechanical stabilities directly from second 
order variations of free energy (J13j) . In the classical literature on stabilities of shells, such 
as Ref.fl and Q, the critical pressure (J2J) are obtained under a special assumption of the 
instable mode that the concave part after instability is the mirror image of its initial one. 
Therefore, the present work implies that we can also analysis the mechanical stabilities 
of solid shells directly from second order variations of free energy (0) without the special 
assumption in conventional literature. 

It is well known that cell membranes contains lipid bilayers and membrane skeleton. Our 
theory of polymer membranes may be applicable for the membrane skeleton because it is also 
a cross-linking structure at molecular levels. In the future, we can turn to the elasticity of 
cell membranes after we fully studied the elasticity of lipid bilayer and membrane skeleton. 
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APPENDIX A: SURFACE VARIATION THEORY 

Here we briefly retrospect the surface variation theory originated in Ref. 3| 



and fully 



developed in Ref. [18j. 

We use a smooth and closed surface M. in 3-dimensional Euclid space E 3 to repre- 
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sent a membrane. As shown in Fig. Q we can construct a right-hand orthonormal sys- 
tem {ei,e2,e3} at any point r in the surface and call {r;ei, 62,63} a moving frame. The 
differential of the frame is denoted by 

dr = wid + u 2 e 2 , 

( A1 ) 

aej = uJij&j (i = 1, 2, 3), 

where u\, uj 2 and U{j = —ojji = 1,2,3) are 1-forms. The structure equations of the 
surface are 

gL>i = co>i 2 Aco> 2 ; (A2) 

du>2 = W21 Awi; (A3) 

uj 13 = auo 1 + buj 2 , = bui + ccu 2 ; (A4) 

duij = u ik Au kj = 1,2,3). (A5) 

n 

Readers should notice that the operator "d" is an exterior differential operator |l^J in this 
paper. The area element, mean curvature and Gaussian curvature are respectively expressed 
as: 

dA = ui Alu 2 , (A6) 
H = (a + c)/2, (A7) 
K = ac-b 2 . (A8) 

If A4 undergoes an infinitesimal deformation such that every point r in M. has a dis- 
placement 5r, we obtain a new surface M! = {r'\r' = r + 5r}. 5r is called the variation of 
surface M. and can be expressed as 

5r = 5it + 5 2 r + S 3 r, (A9) 
6iT = niei (« = 1,2,3), (A10) 

where the repeated subindexes do not represent Einstein summation. Due to the deformation 
of M., ei, e 2 , e 3 also change. We denote the change as 

= Qujej, Vtuj = -O yi . (AH) 

Using the commutativity between Siji^= 1,2,3) and d, we obtain the fundamentally 
variational identities of the move frame 

S\UJi = dfli — 0>2^121> (A12) 
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5iU) 2 = Q1UJ12 — W\£l\V2, 


(A13) 


^113 = CiQi, ^123 — b£li; 


(A14) 


O2UJ1 = ^2*^21 — ^2^221; 


(A15) 


5 2 ^2 = dVt 2 — C<Jif2212, 


(A16) 


^213 = bVl2, ^223 = C^2] 


(A17) 


O3LO1 = ^31^31 — 6^2^321 , 


(A18) 


S3LJ2 — ^3^32 — ^1^312, 


(A19) 


dfl 3 = 1^313^1 + ^323^2 5 


(A20) 


SiuJij = dVLuj + QukWkj — Uik^likj- 


(A21) 
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FIG. 1: Smooth and orientable surface M. we can construct a right-hand orthonormal system 
{ei,e2,e3} at any point r in the surface and call {r; ei, e2, 63} a moving frame. 



TABLE I: The correspondence principle between polymer statistics and path integral method in 
quantum mechanics. 

Quantum Mechanics Polymer Statistics 



time t 

i/h 

Lagrangian L = (m/2)(dr / 'dr) 2 — V[r(r)] 

mass m 
potential V(r(r)) 
propagator K = J exp[(i/ft) J Q * LdT]D[r(r)] 
Hamiltonian H = -|-V 2 + V 
ihdK/dt = HK, (t > 0) 



the number of segments N 
-P = -l/T 

energy E = (r]/2)(dTL n /dn) 2 + U(R n ) 

V = 3/G0&8) 
-U(R n ) 

partition function Z = J exp[— (3 J Edn]D[R n ] 



P 



1 T72 



v z -u 



(l//3)dZ/dN = PZ 
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